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One of the major challenges of nuclear physics is to enlarge the present limits of the 
chart of the nuclides [1, 2]. The experimental progress during last decade has allowed 
to reach the proton drip-line up to charge 80 [3, 4]. In contrast, the situation is quite 
different on the neutron-rich side. Recent studies of neutron-rich light nuclei have shown 
very exciting issues, that new magic numbers might appear and some disappear when 
moving to the neutron drip- line [5]. Near the neutron drip- line, the neutron-matter 
distribution becomes very diffuse and of large size giving rise to "neutron halos" and 
"neutron skins" . These neutron-rich objects have sparked renewed interest in the nuclear 
structure theory [6, 7]. 

In this Letter we present our preliminary results in searching for highly neutron- 
excessive stable light nuclei within Hartree-Fock-Bogoliubov (HFB) approach with 
Skyrme effective interaction, which has the following form 

Vij = toil + xoPa)S{r) + (l/2)ti(l + xiP^)[k'^5{r) + 
5{r)k^] + t2{l + X2Pa)k'5{r)k + {l/6)ts{l + x^P,) 
p"(R)5(r) + iWo[k! X 5(r)k]((7, + Oj) 

where r = - r^-, R = (r^ + rj)l2 , k = -i(vl - V^)/2, k' = - %)/2, 

Po- = (1 -I- (7j(7j)/2. Parameters are given in Table 1. 

We have used the set of parameters Ska [8] and compared the results with the most 



Table 1. Parameters of the Skyrme forces 



Force to ti 


t2 


^3 


Xo Xi X2 Xs Wo 


a 


McV fm-'' McV fm-'^ 


McV fm'= 


MeV fm3 + 3a 


McV fm^ 




Sly4 -2488.9 486.82 
Ska -1602.8 570.9 


-546.39 
-67.70 


13777.0 
8000.0 


0.834 -0.344 -1.0 1.354 123.0 
-.020 0.0 0.0 -0.286 125.0 


1/6 
1/3 


widely used set Sly4 [9]. 


In [10] 


some of us 


shown that for deformed nuclei ^^Mg 


and 



^^^•^""^Si the most satisfactory description of observed spectra comes with the set Ska. 
Pairing effects were included in the standard way with the pairing constant G = 19 /A 
both for protons and neutrons and restricted to the space of bounded one-particle states. 



We carried out the investigations of isotopes ^~^^He, i4-44q ss-soq^ 
analyzed how results depend on forces we have used. For isotopes ^~^^He our results 
matched the known ones [11]. For Helium the last stable isotope with respect to two- 
neutron emission is ^He. Comparison with experiment (see the figures) shows that both 
Ska and Sly4 equally good describe the known experimental data. In our calculations 
with forces Ska we have found a nucleon stable isotope (See Fig. 1). Wc did not 
plot the two-neutron separation energy for because all its neighboring isotopes 
are unstable. With forces Sly4 this isotope appears to be unstable with respect to 
one neutron emission, though the last filled level is close to zero and one can talk 
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Figure 1. Calculated separation energies of one (Sn) and two {S2n) neutrons for 
isotopes 14-44Q -^yj^jj different choices of Skyrme forces compared to the experimental 
data [12]. We did not plot the two-neutron separation energy for because all the 
neighboring isotopes are unstable. 

about "quasistability" in this case. From nucleus to nucleus the situation repeats itself, 
whenever the nucleus is "quasistable" with interactions Ska, then it is "quasistable" 
with interactions Sly4. In all investigated cases the last filled level for nuclei close to 
nucleon stability borderline always had a negative parity. And under "quasistable" we 
mean that this nucleon has a non-zero orbital momentum and the resulting centrifugal 
barrier prevents the neutron from emission at its low energies. 

The obtained data for is given in Table 2. One can see that the values of proton 
and neutron deformation for are negligibly small. 

Table 2. Calculated values of binding energy E, neutron and proton separation 
energy Sp^n, root mean square radii rp,„ quadrupole moments Qp^n and deformation 
parameters /32'" for the stable isotope as calculated with Ska forces. 



E 




Sp 


rn 




Qn 


Qp 






MeV 


MeV 


MeV 


fm 


fm 


e^ fm^ 


e^ fm^ 






168.274 


0.593 


36.822 


4.202 


2.943 


0.031 


0.003 


0.004 


0.004 



The maps of proton and neutron distributions in r, z coordinates (incorporating the 
symmetry of the problem) for and are shown in Fig. 2. One can see although 
the proton "cloud" is expanding it remains coated with the neutron halo which is about 
2 fm thick. 

Our preliminary calculations of nearby isotopes showed, that there are stable 
isotopes around ^°0 among even-even nuclei, namely ^0'42,44j^p with one neutron 
separation energies are respectively Sn = 0.13,0.43,0.1 MeV and ^^'^^Mg with one 
neutron separation energies Sn = 0.8,0.67 MeV (See Figure 3 for the details). These 
stable isotopes were also found with Ska forces and except ''^Mg they lie beyond the 
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Figure 2. The map of proton pp and neutron pn distributions calculated for ""^O (a,b) 
and for ^"O (c,d). The proton "cloud" is expanding with the increasing of the neutron 
number, but it remains coated with the neutron halo with 2 fm thickness. 
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Figure 3. The part of the neutron drip-line. For each value of Z the heaviest stable 
isotope has been presented for experimental data (Audi 2003 [12]) and Sly4 forces. 
For the very neutron-rich isotopes HFB calculations with Ska forces shows the different 
behaviour of the neutron drip-line. See Fig. 1 for the details. 
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conventional stability valley. We have observed that the usual shell closure at N=50 is 
practically absent for the neutron rich Ca isotopes, and at N=40 a new closure appears. 

It should be mentioned, that the stability with respect to neutron emission is defined 
within very narrow range of the binding energy (about 0.5 MeV). Therefore it would be 
desirable to check this effect with other theoretical models. 

Our detailed investigations of the decay properties of the neutron drip-line nuclei 
in this region of possible stabihty (isotopes of C,0,N,Mg with the very large neutron 
excess) are in progress. The results will be published elsewhere. 
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